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We consider the decay of positronium to a neutrino-antineutrino accompanied by a single photon. 

Since the neutrino pair go undetected, this appears as a single photon decay of positronium. These 
decay channel are mediated through the exchange of the massive W and Z vector bosons of the 
electroweak interaction. After summing over the various neutrino channels, the standard model 
calculation yields the rate for such a single photon decay process of T ps^-y = 1.72 xl0“^® s“^. 


I. INTRODUCTION 

Positronium (Ps) is a metastable leptonic atom [T] 
composed of an electron and its antiparticle, the positron. 
This exotic atom resembles in many ways the hydrogen 
atom, but with bound state energy levels appearing at 
half of the energies due to its reduced mass being half of 
the electron mass. The existence of Ps was firstly pos¬ 
tulated by Mohorovicic in 1934 [2] and finally observed 
in 1951 by Deutsch [3]. As in the case of the hydrogen 
atom, two Ps can form a molecular bound state, leading 
to a hydrogen-like molecule but with a leptonic nature 
Ps 2 , which was predicted in the 40’s by Wheeler [4], and 
recently observed Different decay channels of both 
positronium and molecular positronium Ps 2 have been 
addressed in the framework of QED [6H9]. 

Just as the case with the Hydrogen atom, Ps possesses 
two different spin states, a singlet state which is called 
parapositronium (p-Ps), and the triplet state called or- 
thopositronium (o-Ps). Within the strict confines quan¬ 
tum electrodynamics (QED), which preserves charge con¬ 
jugation symmetry, the p-Ps state can only decay into 
an even number of photons, while o-Ps decays to an odd 
number of photons m- The vast majority of the decay 
modes of Ps which are allowed within pure QED, both 
the dominant modes as well as more exotic decay chan¬ 
nels, have been computed [3. Eor the dominant modes 
the calculations have been extended to include higher 
order QED radiative corrections. The inclusion of the 
electroweak interaction which are mediated by the ex¬ 
change of the massive vector W and Z bosons allows for 
additional decay modes such as those involving neutrinos 
which are absent in pure QED analysis. For instance, o- 
Ps can decay into two neutrinos leading to an invisible 
channel [3 n] which has been experimentally investi¬ 
gated [12] ■ However, the influence of the electroweak 
interaction into the decay modes of p-Ps remain unex¬ 
plored. 

In this work, we report on the calculation of the decay 
mode of p-Ps into a neutrino-antineutrino plus a single 
photon. Since the neutrino pair goes undetected, this ap¬ 
pears experimentally as a single photon decay channel. 
The studied decay mode involves the exchange of mas¬ 
sive vector bosons of the electro-weak interaction and as 
such cannot arise within pure QED. The energy of the 


photon emitted in this mode exhibits a continuous spec¬ 
tra of energy between zero and mg. The process can shed 
some additional light on the influence of the electroweak 
interaction in leptonic bound states, as well constituting 
another novel decay channel for the completeness of the 
decay chart of Ps. 

II. SINGLE PHOTON DECAY MODE OF PS 

The annihilation of Ps into a neutrino-antineutrino 
plus a single photon is a consequence of the electroweak 
interaction. Through the exchange of the massive W and 
Z vector bosons, the electrons and positrons comprising 
the positronium can produce a neutrino-antineutrino pair 
which can be accompanied by the emission of the pho¬ 
ton off the incident electron or positron. To secure the 
single photon total rate (here denoted by Ps —>■ 7 ), we 
compute the contributions from each of the possible fi¬ 
nal states involving each of the three species of neutrinos 
(iZe, r'/j, and add the rates incoherently. For the fi¬ 
nal state containing the photon plus the i>e ~ Ve pair, 
the relevant Feynman diagrams with the final state pho¬ 
ton emitted off the incident electron leg are displayed 
in Fig. 1. In addition, there are the analogous graphs 
with the outgoing photon emitted from the positron line 
which are not displayed. The process with the I’eh'el fi¬ 
nal state includes both t-channel W boson exchange (Fig. 
la) and s-channel Z boson exchange (Fig. lb). In ad¬ 
dition, the photon can also be emitted off the internal 

charged W line as shown in Fig. 2, but its contribu- 
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tion to the amplitude is further suppressed by , with 
me and M\y the electron and W masses respectively, 
and is thus neglected. For the and final 

states, only the s-channel Z boson exchange displayed 
in Fig lb contributes. Here again, there are the anal¬ 
ogous graphs with the photon emitted off the positron 
line which are not explicitly displayed. The W exchange 
diagram of Fig. la does not contribute in either of these 
cases. The energy-momentum of the incoming electron 
and positron are labeled as pi , P 2 respectively, while the 
outgoing neutrino and antineutrino are labeled as ps and 
P 4 respectively and the emitted photon 4-momentum is 
denoted by k. Here, the energy-momentum vectors are 
represented as (E,p) and a ~ 1/137 is the fine structure 
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constant. We work throughout using natural units {h = c 
= 1). Throughout the calculation, the tiny mass of the 
neutrinos is neglected and they are treated as massless. 

The decay rate of a bound positronium state into 
each particular mode, ; I = e,v,T, can 

be expressed in terms of the cross section, cr(e+e“ —>■ 
^ = e,/r,T, associated with the process e^e“ —)■ 
viDi"/ ; Z = e,/i,r in which the e“'",e“ constituents of 
the bound state as a free particles. So doing, the decay 
rate for the single photon decay mode of Ps is then given 
by 


Pps—f7 — Tps— I'l'y 

~ ^ ^ 4'U^e/|^Ps(0)| (1) 

where |^Ps( 0 )p denotes the probability of finding the 
electron and positron in the same point of space and 
Vrei is the relative velocity of the colliding e^e“ pairs. 
Here the electron-positron annihilation cross sections 

e~ ^S 




^ 1 f <PP3 

AmlVrel J 2(27r)3|p3| 

d^P4 f d^k 
2(27r)3|^4| J 2(27r)3|fc| 

{2Tr)'^S^*'^ (pi + P2 - P3 - P4 - fc) < |Tp >, 

( 2 ) 

where < |rp >= j is the totally spin 

averaged square of the transition matrix element T ob¬ 
tained by averaging over the initial e~ and e’*' spins (si 
and S 2 respectively) and summing over the undetected 
final state neutrino and antineutrino spins (S 3 and S 4 re¬ 
spectively). Eq. ([^ includes the Lorentz invariant phase- 
space factors, as well as the energy-momentum conserva¬ 
tion (which is guaranteed by the S functions). Ps has a 
binding energy of 6.8 eV, which gives an upper bound 
to the momentum of the leptons involved in the decay 
of Ps. These values are negligible in comparison with 
the rest mass energy of those leptons, m^. Thus, the ini¬ 
tial energy-momentum vectors are taken as (mg, 0 , 0 , 0 ). 
Using this approximation and applying the momentum 
conservation through the delta function, the integral over 
d^P 4 in Eq. ([^ is performed, yielding 


FIG. 1. Feynman diagrams contributions to the annihilation 
process e^e“ —>■ viVi'y, with I — e,p.,T. Both figures (a) 
and (b) contribute to the case I = e, while only figure (b) 
contributes to the cases I = p,T. In addition to these graphs, 
there are also the graphs where the outgoing photon is emitted 
off the positron leg. 



FIG. 2. Feynman diagram contributing to ^(e^e” —>■ r'et'e7 
process which involves the exchange of two W bosons. The 
contribution from this graph is thus highly suppressed and is 
neglected in our analysis. 


1 f d^p3 1 f d^k 
J \P3\ \^3 + k\ J ^ 

S (2me - IpsI - |P3 + k\- \k\^ < ITp > . 

( 3 ) 

This expression can be further simplified integrating over 
the angular degrees of freedom of jpsl by virtue of the 


kinematic constrains present in the delta function result¬ 
ing in 


*^ 6 + 6“—>7 


1 

4m2urei2'^7r4 


Jo W 



d\P3\ < irp >. 


( 4 ) 
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Finally, the angular degrees of freedom of k are inte¬ 
grated, leading to the final form for the cross section for 
electron-positron annihilation into one single photon 


which in the present case represents the probability to 
have a collision per unit of area and energy of the emit¬ 
ted photon. It”s integral over the photon momentum 
magnitude gives the total cross section as 


1 


^e+e“—>-7 




d\k\ 


e-|fc| 


d\P3\ < irp >. 

(5) 


The relevant T matrix element is secured from the 
Feynman graphs of Fig. 1 and the analogous graphs with 
the photon emitted off the positron leg. Since the mo¬ 
mentum transfer is far smaller than the W and Z masses, 
it can be neglected in the W and Z boson propagators. 
So doing and working in Feynman gauge, one finds 


^e+e 



d\k\ 


da{e 


+ e- 


7) 


d\k\ 


( 8 ) 


The differential cross section is displayed in Fig. 3. The 
cross section is relatively flat as a function of the photon 
energy peaking when the energy of the photon is 2me/3, 
which is represented as the red-dashed line in Fig. 3. 
In other words, the most probable energy of the photon 
emitted will be 2me/3. 
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e^Gp 
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( 6 ) 



FIG. 3. Differential cross section for the process e''"e“ —>■ 
7 as a fnnction of the energy of the emitted photon. The 
red-dashed line represents the photon energy leading to the 
maximum probability of collision. 


where the slash notation, i.e., = 'y’^p^, has been em¬ 

ployed and we have summed over the polarizations of 
the final state photon. Here Gp = 

Fermi constant which has the numerical value 1.17x10“® 
GeV“^ and sin^vv — 0.24 is the weak mixing Weinberg 
angle. We employ the Dirac 7 matrix conventions as is 
defined in Ref. m- The diagrams with the outgoing 
photon emitted from the positron leg has the analogous 
structure except for the replacement the fermionic propa¬ 
gator with the propagator . The 

coefficient (5i,e guarantees that the W exchange graph 
only contributes to the Ve^el final state. 

With the amplitude 'T in hand, the spin averaged tran¬ 
sition probability can be extracted by standard tech¬ 
niques. To do so, we employ Mathematica [T3]. However, 
before calculating the total cross section, let us first con¬ 
sider the differential cross section 

dcr(e+e~ -» 7 ) _ dg(e+e~ -» vivg) , . 


Finally, integrating the differential cross section from 
[ 0 ,me] as instructed by Eq.([^ yields the total cross sec¬ 
tion of the decay mode e-|-e“ —> 7 which is then employed 
in the calculation of the decay rate of Ps —7 yielding 


Fp,^^ = |^p3(0)|2 


370.8 X Glmle^ 
' 2 V 


(9) 


Using |'I'ps(0)p = TOga®/( 87 r) gives a decay rate of 
rps _>.7 = 1.72 X 10“^® s“^. This value is roughly the 
same order as for the pp (invisible) mode for o-Ps [3]. 


III. CONCLUSIONS 

Due to the presence of the electroweak interactions, 
additional decay modes become available for the positro- 
nium atom. In the present work, we addressed what is 
effectively the single photon decay mode of parapositro- 
nium, which emerges from the pipi^ ; I = e,p,T final 
states and the non observation of the outgoing neutrinos 
and antineutrinos. These processes become operational 
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due to the exchange of W and Z massive vector bosons. 
The study of this decay aids in elucidating the role of 
the electroweak interaction in the decay of Ps. Since this 
mode involves the emission of a single photon, it has a 
somewhat unique signature which in principle may be 
detected. The calculated decay rate is rps_>..y = 1.72 x 
10 “^® s“^, which is experimentally challenging, is com¬ 
parable to the previously reported electroweak mediated 
decays in orthopositronium but without the photon emis¬ 
sion. 

The decay mode presented will be also present in the 
decay of other exotic atoms, such as true muonium; the 


bound state of a muon and anti-muon. For this bound 
state, the analogous single decay rate for will be at en¬ 
hanced by a factor of ~ 8 x 10® due to the heavier mass 
of the muon with respect to the electron. 
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